Phase transitions from preheating in gauge theories 
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We show by studying the Abelian Higgs model with numerical lattice simulations that non-thermal 
phase transitions arising out of preheating after inflation are possible in gauge-Higgs models under 
rather general circumstances. This may lead to the formation of gauged topological defects and, if 
the scale at which inflation ends is low enough, to electroweak baryogenesis after preheating. 
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Over the past few years there has been somewhat of a 
revolution in our understanding of the dynamics of the 
end of inflation. The traditional picture of reheating aris- 
ing out of the perturbative decay of the inflaton field as 
it oscillates about the minima of its potential has been 
replaced by the possibility of an explosive particle pro- 
duction during an earlier period, known as preheating. 
During preheating, parametric resonance of the inflaton 
field generates very large fluctuations of the scalar fields 
coupled to the inflaton, leading to the production of large 
numbers of particles Q . Due to the weakness of the in- 
teractions, the short-wavelength modes do not thermal- 
ize, and the effective temperature of the long-wavelength 
modes is much higher than in the standard reheating sce- 
nario. This may lead to symmetry restoration and, when 
the Universe cools down as it expands further, a subse- 
quent non-thermal phase transition . The fact that the 
fluctuations produced during preheating have large occu- 
pation numbers implies that they can be considered as 
interacting classical waves, an important result because 
it means that the dynamics of fluctuations during and 
after preheating can be studied using lattice simulations 
H . A concrete example of a non-thermal phase transi- 
tion occuring after preheating was presented in Ref. ||] 
(see also Ref. ||). The phase transition that they found 
is first-order, and depending on the field content of the 
model being investigated topological defects may form, 
an intriguing result as it opens up the possibility that 
inflation can create a defect problem if for example they 
produce gauged monopoles or domain walls [|| || . 

Non-thermal phase transitions may even solve the old 
puzzle of baryon asymmetry in the Universe Al- 
though the baryon number is conserved perturbatively 
in the Standard Model, there are non-perturbative in- 
teractions that violate this conservation law. The rate 
of baryon number violation is extremely low at low ener- 
gies, but it becomes much higher in the high-temperature 
phase of the electroweak theory. Thus it is possible to 
generate the observed baryon asymmetry if, for some rea- 
son, the fields are out of equilibrium at the electroweak 
scale and thermalize to a temperature below T c . This 
could be the case even in the standard big bang cosmol- 
ogy, if the electroweak phase transition were strongly first 
order, but at least in the minimal Standard Model it is 



not, as lattice simulations have shown flQ(| . However, in a 
non-thermal phase transition, the fields are driven out of 
equilibrium by the oscillations of the inflaton, and baryo- 
genesis may be possible, if the reheating temperature is 
much lower than T c 

Despite the exciting possibility of electroweak baryoge- 
nesis, most of the numerical work on non-thermal phase 
transitions so far has concentrated on scalar fields jjl 
or has been restricted to one spatial dimension Ex|. In 
this letter we present results from simulations of the 
Abelian Higgs model in two rather different cases. The 
first case is a direct analogue of the simulations in Ref. [EJ . 
The inflaton itself is charged under a gauge group and 
eventually breaks the gauge symmetry in a first-order 
phase transition. (For simplicity, we use the terminology 
of spontaneous symmetry breakdown, although the gauge 
symmetry is not actually broken in the Higgs phase.) 
The second case is more relevant for electroweak baryo- 
genesis. We show that even with a Higgs mass that is 
compatible with experimental bounds, the transition is 
sharp, and electroweak baryogenesis is therefore possi- 
ble. Although we restrict ourselves to the Abelian case 
in our simulations, we expect that our conclusions apply 
to non-Abelian theories as well. 

The Lagrangian of our model is 



L = —F^F^ + (£>'- 



>D^-\M 2 ~v 2 ) 2 . (1) 



Here the gauge covariant derivative is D^c, 



ieApcf), and F^ = A v -^ — A^ v . The couplings A and 
e are assumed to be small, and we will use A ~ e 2 in our 
estimates. 

Ideally, we would like to study the quantum field the- 
ory defined by Eq. (Q) , but solving for the time evolution 
of even a simple quantum system is a formidable task. 
Therefore we have to resort to the classical approxima- 
tion, which is expected to work as long as the dynamics 
is determined by modes with a macroscopic occupation 
number H. For studying the dynamics, it is convenient 
to fix the temporal gauge Aq = and use the confor- 
mal time rj defined by drj = dt/a and the rescaled fields 
<j> = a<f>, Ei = —dr/Ai. The equations of motion for <f> and 
Ai follow from Eq. (fil): 
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4>+{2\v z a z + d 2 a/a)4>- 2\\<j>\ J 4>, (2) 
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d v Ei = dj Fy + 2eha4>* Difa (3) 
diEi = 2elm0*<9, ) 0. (4) 

The initial conditions for the fields are those produced 
by inflation: the gauge field is in vacuum and the co- 
variant derivatives of the Higgs field vanish. This allows 
us to fix the remaining gauge degree freedom by setting 
initially Ai = and <f> = <f) Q =constant. 

We separate <fi into the homogeneous zero mode <f> and 
the inhomogeneous fluctuations 6<f> = <j>— <j>. We take the 
quantum nature of the system into account by introduc- 
ing small fluctuations for the fields Ai and 6<j> and for 
their canonical momenta Ei and St: = d v 5<f>. The width 
of these classical fluctuations is chosen to be equal to 
the width of quantum fluctuations in the vacuum calcu- 
lated for free fields. We allow fluctuations in the phase 
of 8<p and fix the associated gauge degree of freedom by 
choosing diAi — 0. The longitudinal component of Ei is 
determined from the Gauss law (Q). 

In the very beginning of our simulation, when the fields 
are in vacuum, the conditions required by the classical 
approximation are not satisfied, but we expect that the 
final results will be unaffected. What is important is not 
the precise nature of the initial fluctuations, but that 
some small fluctuations are present. 

As in the scalar theory 0], the time evolution begins 
with a period of parametric resonance. The resonance 
parameter q is given by q « e 2 /A. Let us first consider 
the case q 3> I, in which the resonance is broad, and 
during the first oscillations, a large amount of energy is 
transferred from the zero mode <f> to the long-wavelength 
modes p ~ A 1 / 2 ^ of Ai and S<f>, from which it soon 
spreads to all modes with p<p* ~ e<f> . We can ap- 
proximate the state of the system after this period by 
assuming that the modes with p < p* thermalize to some 
effective temperature T e g, but those with p>p* remain 
in vacuum. Then the energy density in these fluctuations 
is 

„ f* d*p 2 T cS 3 

P ~ ] J2^ P -pr~ p * TcS > (5) 

and after preheating it is of the same order as the ini- 
tial energy density in the zero mode po ~ e 2 cj) , which 
implies T c g ~ 4>ol e - I n the reheating picture, the tem- 
perature after the equilibration of the fields would be 
T r ~ yfetpQ < T off . 

Since the occupation number of the long-wavelength 
modes is n p ~ T c s/p, which is large whenp <p* provided 
that eCl, the classical approximation works well after 
preheating begins. 

The zero mode <f> continues oscillating around the min- 
imum, but the fluctuations in 5(f) and Aj induce an effec- 
tive mass term 

m c 2 ff « -2Xv 2 + 4A<<50 2 ) + e 2 (A 2 ) . (6) 

The magnitude of the fluctuation terms is 



W 2 > ~ (A 2 ) ~ /"* ^r off (7) 

J (27r) J p z 

In the reheating picture, the fluctuation terms would be 
much smaller, (5(f> 2 ) ~ (A 2 ) ~ T 2 ~ ec/) . This shows 
that m 2 s can become positive, thereby restoring the sym- 
metry, even if the reheating temperature is below T c . 

When the Universe expands further, the fluctuation 
terms decrease and the system undergoes a phase tran- 
sition to the broken phase. The nature of this transition 
can be studied by calculating the effective potential of 
<j> in the background of the fluctuations 5<j) and A^. If 
e 2 3> A, the contribution from Ai will be more impor- 
tant. Taking the one- loop contribution from the gauge 
field into account, we have 

V eS $) « -2Xv 2 f + A0 4 + T cS ** log P —^, 

J (2^ p 2 

(8) 

where tua ~ ecf) is the photon mass generated by the 
Higgs mechanism. 

To understand the shape of the potential (||), we ex- 
pand it both for small and large (f>, 

V (6)~f m ° ff ^ ~ Cie - TcS ^ + ? < p */ e ) 

off W ~ | C2Tcgp 3 2Xv 2f + ^ ^ >> ^ /e) 

(9) 

where C\ and C 2 are numerical factors and m 2 ff is given 
byEq. (|. 

The origin <f> = is a local minimum whenever m 2 s is 
positive. Assuming first that p* > ev, the cubic term in 
Eq. ([)]) induces another minimum for the potential when 
m 2 s becomes small enough, and eventually when this 
new minimum becomes the global one the system enters 
the Higgs phase in a first order phase transition. While 
this phenomenon is present also in equilibrium Jl3[ ] , the 
transition is stronger in our case, since the cubic term 
is proportional to T c g ^> T r . The existence of this min- 
imum requires e 2 > A, since otherwise the contribution 
from the scalar loop, which does not contain any cubic 
term, would dominate in Eq. @. 

Even if e 2 < A, the potential can have two minima, pro- 
vided that < ev Then the tree- level minimum 
is the global one if the logarithmic term in Eq. @ is 
smaller than Xv A , i.e. T c gpl<Xv A . It is difficult to si- 
multaneously satisfy this inequality, along with the con- 
dition m 2 s > 0, and we were unable to do this in our 
simulations. 

To confirm this expected behaviour, we carried out a 
numerical lattice simulation with our model. We chose 
~4> a = 0.25Mpi and A = 2 ■ lCT 13 , e = 6.4 • 1(T 6 , 
v = 7.2 ■ 10 _4 Mpi. The lattice spacing was Sx = 
9.3 • 10 5 A/ P1 \ the time step 5f] = 1.2 ■ 10 5 M P " 1 1 and the 
size of the lattice was 320 3 . The universe was assumed 
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FIG. 1. The time evolution of \<f>\ in the first simulation 
on a 320 3 lattice with the initial condition (j> — 0.25Mpi. 
That \(f>\ 2 is below v 2 indicates that the symmetry is restored. 
At r] w 3 • IC^MjTj 1 , the transition to the Higgs phase takes 
place. The inset shows the effective temperature of different 
Fourier modes of the electric field Ei measured at various 
values of r\ (given in units of IC^Mp! 1 ). The energy density 
in the short-wavelength modes is suppressed by a factor of 
10 4 relative to the long-wavelength modes even at the end 
of the simulation, and therefore the discretization errors are 
expected to be small. 

to be radiation dominated with a{rj) = 1 + r]H, where 
H = 8.3-10 _8 Mpi. Since (j> is not a gauge-invariant quan- 
tity and can therefore only be defined in the vacuum, we 
did not measure its value. Instead, we show \(j>\ 2 as a 
function of time in Fig. [j] The fact that |</>| 2 < v 2 when 
l^-l^Mp/ < r) < 5-l0 9 M pi 1 clearly shows that the gauge 
symmetry is restored and the system is in the Coulomb 
phase. The amplitude of the oscillations remains quite 
large, which is probably a finite-size effect. In an infinite 
system, there would be more infrared modes to which the 
zero mode of </> could decay. Eventually, the system en- 
ters the Higgs phase in a first-order transition, as in the 
scalar theory. The first-order nature of the transition can 
be seen from the configurations during the transition; for 
example by looking at the isosurface of \(j)\ 2 we would see 
a growing bubble of the Higgs phase characterized by a 
larger value of \4>\ 2 . 

In order to check that the separation of scales below 
and above p» indeed takes place, we measured the effec- 
tive temperature of different Fourier modes of the elec- 
tric fields Ei at various times during the simulation. A 
reason for choosing this quantity rather than the power 
spectrum of </> or Ai is that it is gauge- invariant. In equi- 
librium |i?i,/t| 2 is constant and its magnitude is propor- 
tional to the temperature. Therefore we can use it to 
define the effective temperature of a single mode 

1 H 3 k 
T off (p = fc/ a ) = -|£ T fe | 2 — -, (10) 
la ' ( 2tt) 
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FIG. 2. The time evolution of \(j>\ 2 in the second simulation 
on a 240 3 lattice with the initial condition (j> = 1 TeV. Again, 
the symmetry is restored, and at rj ~ 3 GeV -1 , the transition 
to the Higgs phase takes place. The effective temperature 
measured at various r\ (given in units of GeV -1 ) in the inset 
shows that in this case the resonance is narrower, but the 
temperature still grows to high values. 

where the superscript T indicates that we have included 
only the transverse component kiEj k ~ 0, since the lon- 
gitudinal component is fixed by the Gauss law. The inset 
of Fig. [I] shows the product aT e g, which is the effective 
temperature of the rescaled fields, as a function of the 
conformal momentum k = pa. 

Immediately after preheating, the temperature of the 
long-wavelength modes is T c g ~ I0 4 Mpi and the occupa- 
tion number n p = T cS (p)/p ~ 10 10 is huge. The cutoff 
momentum is p* — k^/a ss I0 -6 Mpi/ 'a. With time, the 
modes with higher and higher k thcrmalize and the tem- 
perature decreases, but since the couplings are small, this 
process is very slow. The modes with i> fc* are strongly 
suppressed even after the phase transition, and there- 
fore we believe that the lattice approximation remains 
reliable even at the end of the simulation. Because the 
modes with the highest momenta do not remain exactly 
in the vacuum, discretization errors cannot be ruled out 
completely. 

For the electroweak theory, the opposite case q < 1 
is more relevant, since q ~ m^/mjj. In this case, the 
parametric resonance is narrow and the energy transfer 
is less efficient. However, since the expansion rate of the 
Universe is much slower in this case, it could still lead to 
a similar phenomenon. Most of the energy of the infla- 
ton is transferred to a narrow momentum range of the 
fluctuations, but the long-wavelength modes thermalize 
and the energy is spread to all long-wavelength modes. 
After that, the system should behave as in the case with 
a broad resonance. In equilibrium, the phase transition 
is not of first order if e 2 > A, but as discussed earlier, we 
expect the transition to be stronger in our case. 



The realistic values for the couplings in the electroweak 
theory would be A ~ e 2 ~ 1, but in that case our simu- 
lations are not reliable. With these couplings, the inter- 
actions are important even in the vacuum state, and the 
classical approximation cannot be trusted. Therefore, 
we have used slightly smaller couplings, A = 0.04 and 
e = 0.14, which allow us to use the classical approxima- 
tion. The initial value of the Higgs field was <f> = 1 TeV. 
We also chose v — 246 GeV and a{rj) = 1 + tjH with 
H = 0.7 GeV. The lattice spacing was Sx = 1.4 TeV" 1 , 
time step Sr] = 0.14 TeV" 1 , and the lattice size 240 3 . 

In this case, <j> cannot be the inflaton, because its cou- 
plings are much too strong. However, the homogeneous 
initial condition for <f> may arise from a previous pre- 
heating phase, in which (f> couples to the inflaton with 
a coupling constant that is much smaller than e. Then 
the parametric resonance will transfer a large amount of 
energy to modes of <f> with very long wavelengths. The 
alternative possibility is that quantum fluctuations of <j> 
give it a large spatial average during inflation. 

As in the earlier case for the inflaton, we show \4>\ 2 
and the effective temperature of different modes of E.- L 
in Fig. (0). This time, the energy is transferred into 
a narrow band of gauge held modes. Nevertheless, the 
long-wavelength modes thermalize, and we reach a simi- 
lar situation to that in the first case, in which the long- 
wavelength modes k < 300 GeV have an effective temper- 
ature T e ff « 10 4 GeV, and the symmetry is restored. At 
?y w 3 GeV" 1 , the system undergoes a phase transition 
to the Higgs phase. The transition is not of first order, 
but it is still rather sharp. 

In the electroweak theory, the conservation of baryon 
number would be violated by sphaleron configurations 
with a rate r sp h ~ ce^X^/p* (T^| when the symm etry 
is temporarily restored, and as discussed in Ref. |L5| , 
the oscillations of the Higgs field could create a large 
baryon asymmetry. If the transition to the Higgs phase 
is sharp enough, the baryon number violation ceases in- 
stantaneously, and the produced baryon asymmetry re- 
mains. 

Our simulations show that a gauge-Higgs system ex- 
hibits the same behaviour as the scalar model considered 
in Ref. M. The first case we considered shows that a 
non-thermal phase transition is possible if the inflaton 
is charged under a gauge group. Although we restricted 
ourselves to an Abelian model, we believe that the qual- 
itative features of our results would be the same in non- 
Abelian theories. In many models, this phase transition 
would lead to the formation of cosmic strings or other 
topological defects. 

The second case we considered has the qualitative fea- 
tures of the electroweak theory, and we find that the sym- 
metry gets restored although the parametric resonance 
is narrow, provided that the expansion of the Universe 
is slow enough. Unlike in the standard thermal phase 
transition scenario, the transition to the Higgs phase is 



sharp, which makes it possible to preserve the produced 
baryon asymmetry. This supports the picture of elec- 
troweak baryogenesis at preheating. 
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